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ABSTRACT 

During a scan of the archival BATSE daily records covering the entire 9.1 years (TJD 8369-11690) of the 
BATSE operation, 3906 gamma-ray bursts (GRBs) have been detected. 2068 of these GRBs are previously known 
BATSE triggers while 1838 of them are new non-triggered bursts. It is important that all events were detected 
in the same type of data and were processed with the same procedure. Therefore these 3906 GRBs constitute a 
uniform sample. We have created a publically available electronic data base"^ containing this sample. We describe 
the procedures of the data reduction, the selection of the GRB candidates, and the statistical tests for possible 
non-GRB contaminations. We also describe a novel test burst method used to measure the scan efficiency and the 
information obtained using the test bursts. Our scan decreases the BATSE detection threshold to ^ 0. 1 phs"' cm"^. 
As a first result, we show that the differential logA^-logP distribution corrected for the detection efficiency extends 
to low brightnesses without any indication of a turn-over Any reasonable extrapolation of the new logA^-logf 
to lower brightnesses imply a rate of several thousands of GRBs in the Universe per year. 

Subject headings: gamma-rays; bursts - methods; data analysis 



1. INTRODUCTION 

The sample of gamma-ray bursts (GRBs) detected by the 
Burst and Transient Source Experiment (BATSE) (Fishman 
et al. 1989) onboard the Compton Gamma-Ray Observatory 
(CGRO) is a few times larger than the yield of all other experi- 
ments that have detected GRBs. Despite that the breakthrough 
in 1997 concerning the distance to GRBs was associated with 
the precise localization of a small number of GRBs by Beppo- 
SAX, the value of the large BATSE sample cannot be overesti- 
mated. Problems where a large statistics of GRBs as well as the 
wide brightness range of the BATSE sample are crucial are; 

1 . estimating the total rate of GRBs in the Universe, 

2. searching for different subclasses of GRBs, 

3. observing gravitational lensing of GRBs, 

4. searching for repetition of GRBs, 

5. estimating the intrinsic luminosity function of GRBs. 

The list can be made longer Here, we mention only those prob- 
lems which are directly associated with the goals of this work. 

Nevertheless, the BATSE sample (i.e., the sample of the trig- 
gered events included in the BATSE catalogs'') is smaller than 
what it could be at the actual sensitivity of BATSE. The reason 
for this is a difficult variable background. The trigger thresh- 
old must be higher and the trigger integration time must be 
shorter than for the case of a constant background. Otherwise 
a high rate of triggers would have caused severe problems in 
the data acquisition. The BATSE trigger adjusted to the back- 
ground conditions missed many weak but still highly significant 
GRBs. Some bursts were also missed due to other reasons (data 
readouts, large background, etc.). 

Many of these non-triggered GRBs can be confidently iden- 
tified in the BATSE daily records which cover the whole period 



of the CGRO operation. Off-line searches for non-triggered 
bursts with less rigid off-line trigger criteria can substantially 
increase the sensitivity of BATSE and extend the sample of 
GRBs both in number and in brightness range. 

The first sample of non-triggered GRBs was published by 
Rubin et al. (1992). A systematic search for non-triggered 
GRBs in 6 years of CG7?(9/BATSE data was performed by 
Kommers et al. (1997, 1998,'' 2000, hereafter K97, K98, KOO). 
Schmidt (1999) performed a scan of the BATSE daily records 
for ~6 years (TJD 8365-10528) using a triggering procedure 
close to the BATSE trigger and found about 400 bursts missed 
by BATSE itself. 

The main feature of our off-line scan of the BATSE data is 
the measurement of the efficiency of the GRB detection as a 
function of brightness using artificial test bursts (e.g.. Stern et 
al. 2000b). This method plays the role of being a calibration of 
the experiment. 

In this paper, we concentrate on a description of the data 
reduction and the tests for possible non-GRB contaminations. 
We also describe the publically available data base containing 
our sample of GRBs. Earlier preliminary versions of this data 
base and some first results were briefly reported in Stern et al. 
(1999a, 1999b, 2000a, 2000b). 

In §§2-5, we describe the procedures of the scan, the iden- 
tification of GRBs, and the GRB data reduction. In §6, we de- 
scribe the general characteristics of our GRB sample, compare 
the sample with the BATSE and the Kommers et al. catalogs, 
and describe the data archive containing the sample. Section 
7 is devoted to various statistical tests for possible non-GRB 
contaminations using hardness ratios, angular and latitude dis- 
tributions. We obtain upper limits to the contaminations caused 
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by terrestrial phenomena, solar flares, and Galactic X-ray ob- 
jects. In §8, we estimate the efficiency of the scan and present 
the final logA^-logP distribution, which considerably differs 
from those of the BATSE and the Kommers et al. catalogs. 

2. THE DATA SCAN 

We used the BATSE daily records (DISCLA) from the ftp 
archive at the Goddard Space Flight Center.^ The main array 
of the raw data consists of the number counts in the 8 BATSE 
Large Area Detectors (LADs) in the 4 energy channels with 
1.024 s time resolution for the whole period of observations 
(excluding some data gaps). 

Examples of data fragments are presented in Figure 1 . The 
data are difficult to process because of the large diversity of in- 
terfering phenomena: the variable background, the flaring iono- 
sphere, solar flares, variable astrophysical sources, and occul- 
tations of sources by the Earth causing steps in the count rate 
curves. In addition, the situation is confused by numerous data 
gaps. 
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Fig. 1 . — Two fragments of the BATSE daily records, each corresponding 
to one orbital period of CGRO. Counts rates averaged over ten 1 .024 bins are 
shown for 2 detectors in the 50-300 keV range. Features of the count rate 
curves are: (Upper panel) <51 800 s - strong noise from Cyg X-1 which is in 
its bright hard state; 51 800 s - descent of Cyg X-1; 52 400 s - CGRO reaches 
its most Southern latitude; 54 050 s - rise of Cyg X-1; 54 660 s - particle 
precipitation visible in all 8 detectors; 54 930 s - a solar flare; 55 200 s - 
CGRO reaches its most Northern latitude; 56 270 s - rise of the Crab nebula. 
(Lower panel) 78 500 s - particle precipitation; 79 550 s - descent of Cyg X-1; 
80 250 s - a non-triggered GRB (the brightest and the third brightest detec- 
tors); 81 700 s - rise of Cyg X-1; 83 030 s - a non-triggered GRB (the two 
brightest detectors). This is a rare case of observing two non-triggered GRBs 
during one orbit. Both are very confident. They have peak fluxes of 0.14 and 
0.21 ph s"' cm"^, respectively. 

The complicated background makes usual statistical esti- 
mates of the detection efficiency useless. The estimate derived 
considering only Poisson fluctuations does not work when non- 
Poisson variability of the background is strong. Then there are 
a variety of reasons for missing a detectable burst: a data gap, 
interference with a variable source, or high ionospheric activ- 
ity. The corresponding probabilities are in all cases brightness 

'Available at: http://cossc.gsfc.nasa.gov/pub/data/batse/daily/ 



dependent. In this case, a realistic estimate of the detection ef- 
ficiency should be based on some kind of simulation of the real 
detection conditions using real data. Such a simulation was 
implemented in the form of a test burst method (Stern et al. 
2000b). 

Artificial test bursts were prepared from a sample of 500 real 
triggered BATSE bursts of durations longer than 1 s and were 
added to the original "daily" count rate records. Each test burst 
was created by sampling one of the 500 bursts with a random 
number and rescaling its amplitude to a randomly sampled ex- 
pected peak count rate with a proper Poisson noise. The loca- 
tions of the test bursts were sampled isotropically over the sky 
above the Earth horizon, and their arrival times were distributed 
randomly with an average interval ^ 25000 s (in total, there 
were 11 100 test bursts). The distribution of peak count rates 
prescribed to the test bursts is shown in Figure 2. The shape of 
the distribution is complicated because initially we did not have 
a clear optimal scheme of brightness sampling of test bursts. 
Initially, we sampled the expected peak count rate between 
0.08 counts s~'cm"^ and 0.3 counts s~'cm~^ (uniformly in log- 
arithmic scale) for 60% of the test bursts. This interval was 
believed to cover the main variation of the efficiency curve and 
no detections below 0.08countss~'cm"^ were expected. The 
remaining 40% of the test bursts were sampled in the bright- 
ness interval 0. 3- 3. Ocountss"'cm"^ (according to a power-law 
dN /dP ex P~' '') in order to check whether there exists any ef- 
ficiency versus brightness dependence for bright GRBs. Later 
we improved the efficiency of the search and found that we de- 
tect some real events below the above threshold. Therefore, 
we reduced the lower limit to 0.06 counts s"' cm"^ and later to 
0.05countss~' cm"^. As a result, we have an artificial and non- 
optimal distribution of test bursts. However, it is still suitable 
for estimating the efficiency. 
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Fig. 2. — Peak count rate distribution of test bursts added to the data 
(for transformation of these units to phs"' cm"^ or to total counts, see § 5.2). 
Test bursts are included here independently of whether they were detected or 
missed. The BATSE threshold is around 0.15countss"' cm"^. The shape of 
the distribution is discussed in the text. 

The procedure of the scan is described in Stern et al. (1999b, 
2000b). Briefly, it consists of adding test bursts to the data, the 
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triggering, an express-analysis, the identification, and the clas- 
sification. 

The off-line trigger consists of three criteria to be satisfied 
simultaneously. The first is a traditional one: a significant rise 
in the counts over the estimated background. For the brightest 
detector, a Aa excess was required, and for the second brightest, 
a 2.5(7 excess. This is lower than the BATSE threshold and sim- 
ilar to that of KOO. The excess was checked for in time intervals 
(i.e., triggering time scales) of 1 bin (L024 s), 2 bins, 4 bins, 
and 8 bins. The longer triggering integration time scales gives 
the main gain in the sensitivity as compared to BATSE. The sec- 
ond criterion is a sufficiently high variability around the trigger 
time. The variability is expressed as the residual x^ after mak- 
ing a linear fit of the signal. This x^ should exceed a threshold 
value, 2.5 per degree of freedom. This value was adjusted em- 
pirically in a trial scan. The third criterion is a check whether 
the signal can be attributed to Cyg X-L The signal was fitted 
with the photon flux from the Cyg X-1 direction and the resid- 
ual was checked for a sufficient variability as in the previous 
case. The two latter criteria, which were not used previously, 
turned out to be very efficient in eliminating false triggers re- 
ducing their number by more than an order of magnitude. 

Each trigger was followed by a preliminary estimate of 
whether the trigger is a GRB candidate (the person performing 
the scan was unaware of whether the candidate event was a real 
or a test burst). Most triggers were rejected by a visual analysis. 
Usually false triggers were caused by occultation steps, sharp 
variations of the ionospheric background, or solar flares. The 
selection criteria used at this step are similar to those described 
in §4, they are just not as stringent, and the count rate was just 
assumed to be proportional to the detector projection area in the 
burst direction (ex cos9). For all selected events, we recorded 
fragments of the original data and some preliminary estimates 
of their parameters. 

At the next step, we identified which recorded events were 
test bursts or BATSE triggers. Finally, we processed all 
recorded events and made a more careful selection of GRBs 
using the procedures described in §3 and §4. 

The rate of our off-line detections of GRBs depends on time 
(see Fig. 3). It reflects the quality of the data, the intensity of 
noise-generating sources (Cyg X-1, X-ray novae), and the so- 
lar activity. One can see a gradual increase in the number of 
detected GRBs with time. It results from improvements in the 
software and from the experience gained during the scan. Note 
that our rate of observing non-triggered GRBs is well corre- 
lated with that of K98. The scan itself took for our group 2 
years. The speed of the scan gradually increased. Now the 
primary scan would take about 200 full working days for one 
person and a comparable amount of work would be spent on the 
second stage data processing. 

3. THE FITTING OF THE GRB LOCATION AND PHOTON FLUX 

The aim of the fit is to find the best burst location, spectrum, 
and intensity describing the count rate data in 8 detectors and 
4 energy channels. This is done using the detector response 
matrix (hereafter DRM), D/jk, where I is the index of the pho- 
ton energy bin, k is the channel number, and j is the detector 
number. D/jk depends on the burst location, the satellite orien- 
tation, and the location of the Earth with respect to the CGRO 
coordinate system. To calculate the DRM we used a version of 
the code written by G. Pendleton which was used for the local- 
ization of the BATSE triggered bursts (Pendleton et al. 1999, 
hereafter P99). The localization procedure itself differs from 



that implemented in the LOCBURST code by P99 in some de- 
tails as described below. 
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Fig. 3. — Time history of the rate of GRB detections during the scan for 
(from top to bottom) test bursts, all real bursts, non-triggered bursts, and non- 
triggered bursts of KOO. The two deepest minima at TJD 8800-8900 and 9900- 
10000 are caused by bright X-ray novae. 

As a spectral hypothesis, we used the Band parametrization 
(Band et al. 1993), which consists of 4 parameters: the low en- 
ergy power law slope, a; the high energy power law slope, f3; 
the transition (peak) energy, Eo; and the normalization factor. 
The problem is that we have only 4 experimental values: the 
counts in the 4 energy channels. Therefore a straightforward fit 
using the Band function would be degenerate. 

We solved this problem using a sample of Band parametriza- 
tions (Band et al. 1993) of real time-integrated spectra of 54 
GRBs. At the first step, we determined which of the 54 spec- 
tra (i.e., which combination of a, j3 and £■(>) that provides the 
best fit of the time-integrated brightest part of the burst using a 
pre-estimated burst location. The location fit used by P99 as- 
sumed that the GRB spectra are power laws (i.e., a 2-parameter 
description), but this is a worse description of real GRB spectra 
as compared to the Band function. 

For a given spectrum Fi we have 



^'7* = ^^DijkFiIi/F^Q-jQQ = djkli, 



(1) 



where C^i is the net GRB signal in counts in the i-th time bin, 
the j-th detector, and the k-th energy channel, F50-300 is the in- 
tegral of the incident photon number spectrum in the 50-300 
keV energy range, and djk is the spectrum-integrated detector 
response matrix. In this normalization, /, is the photon flux in 



units of phs ' cm ^. The value to be minimized is 



ijk 



(djkh + Bijk-Sijk) 

Sijk 



(2) 



where Bijk is the estimated background for time bin / and Sijk is 
the measured count rate. The following parameters are free to 
vary. For the whole event there are three parameters: two pa- 
rameters to describe the burst location (the factor djk depends 
on location) and one parameter for the choice of the spectrum 
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from the Band sample. For each time bin there are two more 
parameters: the photon flux /, and the spectral peak energy Eq 
(at this stage, we kept a and /3 fixed at their best-fit values). 
Having two parameters (/,■ and Eq) free at every time bin is sta- 
tistically meaningful and our chosen model describing the GRB 
spectral evolution as a variation of Eq seems natural. 

The background is estimated as a linear function of time, in- 
dependently for each detector and each energy channel using 
two fitting windows, one before and one after the main peak of 
the event. These windows were set manually using a visual in- 
terpretation of the count rate time profile. In some rare difficult 
cases, only one background fitting window was used. 

The value of x^ was calculated for a few thousand locations 
over a wide area of the sky. Using the x^ map of the sky, we 
defined the best fit location, the 1 cr error boundary, and the qual- 
ity of the event as a GRB candidate (see §4). The most serious 
problem in the location fit is the poor accuracy of the DRM for 
the detectors which do not see the burst directly but still detect 
the burst due to photon scattering by the Earth's atmosphere or 
due to penetration through the satellite. A precise calculation 
of the DRM is an extremely difficult problem as one should 
take into account the 3D mass distribution of the CGRO. In the 
LOCBURST algorithm, the final localization fit is performed 
only for the 4 or 6 brightest detectors. This reduces the system- 
atic error for the best fit location. 

However, this option causes a problem when we make the 
localization for a weak burst where a wide x^ map of the sky is 
needed (see §4). For different locations, different detectors be- 
come the brightest and it is not easy to find the 4 or 6 brightest 
detectors from the data. As far as our main interest concerns 
weak events which requires a large area x^ map, we preferred 
to use the fit for all 8 detectors. Another advantage of the 8 
detector fit is a slightly better signal to noise ratio for the fitting 
net signal given by equation (1). 

The systematic errors of our procedure is probably slightly 
worse than that of LOCBURST. However, the difference is con- 
siderable only for very bright bursts. A more exact localization 
could be achieved, in principle, by iterations using a refined 
DRM and reducing the number of detectors participating in the 
fit. Since the exact localization of strong bursts is beyond the 
scope of this work at its present stage, this was not done. For 
deviations between our and the BATSE locations, see §6. 

The best fit value of the photon flux, /,, for a given burst lo- 
cation and photon spectrum can be obtain by minimizing the 
expression for xj similar to equation (2) with respect to /,, i.e. 
by solving Oxf/dl, = 0: 

_ Y.ik^Sijk-Bijk)djk/Sijk 

J2jkd%/Sijk 

The fitting count rate (that we prefer to use instead of the pho- 
ton flux, /,, for estimating the peak brightness) is then defined 
by equation (1), or, in units of counts s"' cm"^ as 

Cik = (dmk/Sm)Ii, (4) 

where m is the index of the detector having the largest projected 
area Sm in the burst direction. 

The quantities used in this paper to represent the time pro- 
files and the brightness of GRBs are: 

1 . The count rate in the 50-300 ke V energy band and the i-th 
time bin, c,- = cq + ca counts s"' cm"^, with cq and c,3 given by 
equation (4). 

2. The peak count rate, c counts s"' cm"^, which is a function 
of c,, see §5.2. This quantity is also denoted as P (mostly on 



(3) 



figures) following the tradition. 

3. The net (fitting) count rate, C,i = '^jCijk counts s"\ where j 
is the detector number, see equation (1), or the net count rate in 
the 50-300 keV energy band: C, = Cq + C,3. 

4. The peak photon flux, / phs"' cm"^, which is a function of 
/;, see §5.2. 

4. THE IDENTIFICATION OF GRBS 

The main quantitative criterion for selecting real events as- 
sociated with astrophysical sources from the variations of the 
terrestrial background or other fluctuations was based on the 
residual x^ map of the sky. The same criterion was used as the 
significance threshold. Calculating the residual x^ as a func- 
tion of the burst location (eq. 2) we define the "Icr area", i.e., 
a region where x^~Xo < \/2]V where A^ is the number of the 
degrees of freedom, and Xo is the minimum value of x^- The 
location error according to our preference is the maximum dis- 
tance, di, from the best fit location to the boundary of the Icr 
area. Similarly, we define the 4cr region and the maximum dis- 
tance, S4, from the best fit location. Our criterion which any 
event without any exception should pass to be considered as 
a GRB candidate is 64 < 90°. Note, that 4<t is the excess of 
the residual x^ for the fit at a 90° displacement. The statisti- 
cal significance of an event (defined by the x^ with respect to 
the background) is always above 7<t when the S4 criterion is 
satisfied. The 64 criterion efficiently rejects the fluctuations of 
the diffuse ionospheric background. Examples of events which 
have 64 close to its threshold value are shown in Figures 4 and 5. 

Further criteria were more qualitative and subjective. First, 
we checked whether one of the variable sources (the Sun, Cyg 
X- 1 , or one of the known X-ray transients currently active) was 
inside or close to the la area. Then we applied additional re- 
quirements. First, a large hardness ratio (number of counts in 
the 50-100 keV energy range exceeding that in the 30-50 keV 
range), especially for the events close to the position of the 
Sun. Second, for events close to the positions of the hard X- 
ray sources such as Cyg X-1 or X-ray transients, a large dom- 
ination (estimated visually) of the event over the typical cur- 
rent fluctuations of the source was required (a large domination 
over both Poisson and non-Poisson noise was a general require- 
ment). Short soft events were identified as outbursts of X-ray 
pulsars and their locations were always consistent with the lo- 
cation of one of the known X-ray pulsars. 

For each event we checked visually the signal in different 
detectors with the requirement that the signal is visible in more 
than one detector However, no requirements for the signifi- 
cance of the signal in the second brightest detector were applied 
- the signal should just be visually recognized. 

Special attention was paid to the "context" of the event: a 
GRB should be well isolated in order to be accepted. If a persis- 
tent noise from approximately the same direction and of com- 
parable intensity or a high ionospheric activity were seen, the 
event was discarded. The isolation criterion is a subjective one 
to some extent as it is very difficult to formalize. 

There were no "good" events with a localization significantly 
below the Earth's horizon (all tests bursts were sampled above 
the horizon). Such events were only marginally significant or 
did not satisfy the (54-criterion. Therefore there was no reason 
to use the location relative to the Earth's horizon as an inde- 
pendent criterion. Nevertheless events located near the horizon 
were treated more suspiciously if the ionospheric activity was 
high. 
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Fig. 4. — Residual y^ maps of the sky for three weak GRBs (09394b, 
09575a, and 09590c, where the name consists of TJD and an identifying let- 
ter), which passed the selection criteria being near the &i, threshold of the sig- 
nificance (see text). These events were chosen among the least reliable cases. 
The time profiles for these events are shown in Figure 5. Isocontours show 
constant levels of the residual x^ measured as (x^ — Xq)/"" where Xo i^ the 
minimum, a = \/2N, where N is the number of the degrees of freedom. The 
levels at 0.5(T, la (bold contour), 2cr, 3(T, Act (bold contour), 6(T, and Sct are 
shown. The dashed area represents the Earth. The cross shows the location of 
CygX-1. 

5. ESTIMATING THE PEAK BRIGHTNESS OF GRBS 

5. L The Choice of the Measure for the GRB Brightness 

There are two natural choices for the measure of the strength 
of GRBs. One is a time integral of the GRB signal: the photon 
fluence, the energy fluence, or the total counts. Another option 
is to use the peak characteristics: the peak photon flux, the peak 
energy flux, or the peak count rate. If we try to estimate the 
GRB distance distribution, a measure which is closer to a stan- 
dard candle (i.e., that has a smaller intrinsic dispersion) would 
be preferable. We do not know if this is the case - probably 
both the fluence and the peak flux have very large dispersions. 

The choice is, however, natural as the peak characteristics are 
much easier to measure. Besides, the detection efficiency de- 
pends on the peak count rate rather than on the time-integrated 
signal. Among the peak characteristics we choose the peak 
count rate as a measure of the GRB brightness for the following 
reasons. First, the peak count rate unlike the peak photon or en- 
ergy flux does not depend on the assumed photon spectra which 



could be rather uncertain. It is defined by direct measurements. 
Second, the detection efficiency is a direct function of the peak 
count rate, but not of the peak photon flux. The definition of the 
peak count rate when we have 8 detectors with different orien- 
tations. The best choice at this step is Xht fitting count rate, Cik 
or its reduced value c,i (see §3), i.e., the best fit hypothesis of 
the true signal. 
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Fig. 5. — Time profiles for the 3 events that are shown in Figure 4. Fit- 
ting count rates, Ca + Ca (for channels 2-1-3, see section 3), with subtracted 
background are given. 

5.2. The Estimate of the Peak Count Rate 

The weakest GRBs we deal with have a peak count rate of the 
same order of magnitude as the Poisson fluctuations at 1 .024 s 
time resolution. Therefore, if we define the peak count rate as 
the content of the brightest time bin, we probably just choose 
the highest Poisson fluctuation and thus overestimate the bright- 
ness, perhaps by a factor 2 or more for the weakest events. If 
we smooth the signal, e.g. by averaging over several time bins, 
we reduce the Poisson noise but then we may lose true short 
peaks of the GRB signal. A compromise between smoothing 
and preserving the time resolution is the two-iteration scheme 
described in Stern, Poutanen & Svensson (1999). 

In the first iteration, we find the shortest time scale, Ar, = 
2J X 1 .024 s, y = 0, 1 , 2, 3, on which the signal has a significant 
variation between neighboring time bins (using a 4.3a thresh- 
old). Then we search for the Ar,-interval where the count rate 
is the largest. (If AT = 1.024 s, the search is completed.) If 
j > 0, we make a second iteration searching for a significant ex- 
cess in a fraction STi of the brightest Ar, -interval (I < j). The 
significance threshold, h/, depends on the time scale 2' x 1024 s 
at this step. The values of hi were optimized empirically com- 
paring prescribed (expected) and measured count rates for test 
bursts and were set to ho = 3.5<t, hi = 2.4(7, hj, = 1.4(7. 

The result of this method is shown in Figure 6. One can see 
that a reasonable linearity between the expected and the mea- 
sured count rates has been achieved. Nevertheless, after getting 
rid of the Poisson bias we still have a smaller bias of a similar 
origin. Events enhanced by positive Poisson fluctuations have 
a higher probability to be detected. Therefore we observe some 
systematic excess of the measured amplitude over the expected 
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one for the weakest bursts where the detection efficiency is low. 
This bias is accounted for by the efficiency matrix (see §8.2). 
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Fig. 6. — Measured vs expected peak count rates for detected test bursts. 
Several points having a measured count rate much lower than the expected one 
are affected by data gaps causing the peak of the event to be lost. Two points, 
where the measured count rate is much higher, are probably due to confusion 
of the test burst with an overlapping brighter real event. 

It is useful to give the approximate ratio between the peak 
count rate and the peak photon flux. This ratio versus the peak 
count rate is shown in Figure 7. There is no evident correlation 
between the ratio and the brightness. The relation between the 
peak count rate and the peak photon flux can be expressed as 
c = 0.7504 /± 0.063, where the rms variance is given as the er- 
ror. Note that in rare cases the 50-300 keV count rate exceeds 
the 50-300 keV photon flux. The main reason for this is the 
atmospheric scattering of photons. The photon flux by defini- 
tion is the flux of direct photons, while the photons scattered 
in the Earth atmosphere contribute to the count rate. For some 
geometries, the scattered component is considerable. Another 
effect is that photons harder than 300 keV can give a signal in 
the 50-300 ke V range as they interact in the detector via Comp- 
ton scattering. Both effects enhance the 50-300 keV count rate 
with respect to the photon flux in the same energy range and are 
more pronounced for harder GRBs. 

6. THE SAMPLE AND THE DATA ARCHIVE 

We performed the scan of the DISCLA data for the full 9.1 
years of BATSE observations up to the CGRO deorbiting (TJD 
8369-11690). The number of events in our sample classified 
as GRBs is 3906, 1838 of them are non-triggered, and 2068 
we identified with BATSE triggers. The peak count rate distri- 
butions of these GRBs are shown in Figure 8. Note that some 
non-triggered events are very strong (the strongest one has a 
peak flux of 24phs"' cm"^). The reason for BATSE not detect- 
ing some strong bursts is the dead time when the trigger was 
disabled during data readouts or when CGRO passed through 
regions of very high ionospheric activity. On the other hand, a 
few brightest BATSE GRBs are missing from our sample be- 
cause their amplitude overflowed the two byte integer in the 
BATSE archive records and we were not able to process them. 
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Fig. 7. — Ratio of the peak count rate in channels 2+3 (counts s"' cm"^) 
to the peak flux in the 50-300 keV band (phs"' cm"^) versus the peak count 
rate. A few events with c/I > 1 are due to strong atmospheric scattering and 
to the hard part of the spectra contributing to channels 2 and 3 due to Compton 
scattering in the detectors. 

Durations of GRBs were estimated using Tgo (Kouveliotou et 
al. 1993), i.e., as the time interval between the emission of 5% 
and 95% of the total burst photon fluence. The measurement of 
T<)() has considerable uncertainties when applied to weak bursts. 
The result depends on the time intervals (windows) where the 
signal and the background are measured and therefore depends 
on the subjective impressions of the researcher estimating what 
is the signal and what are the fluctuations of the background. 
The duration distribution for our sample is shown in Figure 9. 
One can see that non-triggered bursts are slightly shorter on av- 
erage. This can be explained by the brightness dependent bias: 
losses of some episodes of weak bursts. 
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Fig. 8.— Peak count rate distribution, rfW/rfClogP), of the 2068 BATSE- 
triggered GRBs detected in our scan {thick line histogram) and of all 3906 
GRBs detected in our scan {thin line histogram). The distributions have not 
been corrected for the efficiency. 
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Fig. 9. — Duration distribution (Tgo) for 1965 triggered (thin line histogram) 
and 1713 nontriggered (thick line histogram) bursts detected up to TJD 1 1499. 
The distribution at Tjo < 4 s is biased (see text and Fig. 10). 
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Fig. 10. — Comparison between the BATSE and our estimates of Tjo for the 
1965 bursts detected both by BATSE and by us up to TJD 1 1499. For clarity, 
to avoid the 1 .024 s discreteness of the data, our estimates are dispersed within 
1.024 s bins (otherwise points would merge into horizontal lines). 

The relation between our and the BATSE estimates of Tgo is 
shown in Figure 10. Large deviations are associated with events 
consisting of widely separated episodes with significantly dif- 
ferent brightnesses. One of the two teams can easily lose the 
weakest episode. Another problem is poor precision in the es- 
timation of the fluence distribution as it depends on the details 
of the background fitting. If there is a precursor or an aftercur- 
sor with the fluence around 0.05 of the total one, then one team 
can include this episode into Tgo, the other team can exclude 
it. Note, that deviations between our and the BATSE estimates 
are of both signs and approximately symmetric for long GRBs. 
Our estimates of Tgo for short events are strongly biased due to 
the 1 .024 s resolution and the Poisson fluctuations. There is a 



large probability that an event with Tgo < 1 s will be distributed 
in two time bins (and therefore will be estimated as 2 s long) or 
will merge with a nearby Poisson fluctuation. 

The distribution of the best fit locations of the GRBs in our 
sample over the sky is shown in Figure 1 1 . No significant de- 
viations from the BATSE exposure function (see Paciesas et al. 
1999; and Figure 17) have been found. 
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Fig . 11. — Sky distribution of all 3906 GRBs on an Aitoff-Hammer projec- 
tion in Galactic coordinates. 



6.1. Comparison with the BATSE and Kommers et al. 
Catalogs 

During the period of our scan BATSE triggered 2704 GRBs. 
Therefore we missed 636 of the BATSE triggered GRBs. We 
estimated at a previous stage of the scan (Stern et al. 2000b) 
that ^ 70% of these GRBs were lost due to gaps in the DISCLA 
data, ^ 20% were too short to be detected at 1.024 s resolution, 
and ~ 10% were missed due to various human mistakes. 

K98 scanned the time interval TJD 8600-10800 and found 
873 non-triggered GRBs. We found 1 132 non-triggered bursts 
during that time interval, 745 of them are in the catalog of 
K98. Kommers (1999, private communication; see also KOO) 
inspected the 387 of our non-triggered events that K98 missed. 
Kommers confirmed 224 of these events as probable GRBs. 24 
of the events were classified by him as particle precipitations, 
and 7 as noise. 90 of our events were not classified as the off- 
line trigger of K98 missed those events due to nearby data gaps. 
31 events were classified as "unknown" for their softness (i.e., 
a very weak or no significant signal in energy channel 3). 

We checked most of the 128 events from the K98 catalog that 
are missing in our sample. We confirm 90 of them as GRBs, 13 
we classify as particle precipitations, 4 as long delayed aftercur- 
sors of triggered GRBs, and 2 as solar flares. 12 of these 128 
events were missed by our trigger. We do not include the 90 
GRBs of K98 confirmed by us into our sample as the efficiency 
we estimate does not account for such additional events. 

The comparison of our estimates of the peak fluxes of GRBs 
in the 50-300 keV range with those of K98 and BATSE is 
shown in Figure 12. One can see a reasonable agreement ex- 
cept for a few events where the disagreement could result from 
errors in the setting of the fitting windows. This could, for ex- 
ample, cause part of an event to be confused with the back- 
ground. In general, the agreement between our and the BATSE 
estimates is better than between our and those of K98 in the 
medium brightness range where a meaningful comparison can 
be made. 
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Fig. 12. — Our estimates of the peak flux in the 50-300 keV range versus 
the K98 estimates for the same non-triggered bursts (crosses) and the BATSE 
estimates for triggered bursts (dots) in 1 .024 s time resolution. For clarity, only 
first ~ 30% of the sample is shown. 
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Fig. 13. — Angular deviations between our and the K98 best fit locations 
(crosses) and between our and the BATSE locations (dots) versus the la local- 
ization errors according to our estimates. For clarity, only first ~ 30% of the 
sample is shown. 

The deviations between our best locations and those from 
the K98 and the BATSE catalogs are shown in Figure 13. Most 
of the deviations are within the \a errors. There is, however, 
a substantial tail to much larger deviations caused by system- 
atic errors. Very large deviations can result from the choice 
of different local minima in the location fit. The peak in the 
distribution of the deviations between our and the BATSE lo- 
cations is between 2 and 2.5 degrees for bright GRBs. This 
value approximately characterizes our typical systematic error 
for the locations. The average systematic error of LOCBURST 
as estimated by P99 is ^ 1 .7 degrees. Our procedure is not opti- 
mized for precise localization of strong bursts, since we do not 



take into account nonlinear effects occurring at very high count 
rates and do not turn off less illuminated detectors in the fit as 
LOCBURST does. Our scheme is adjusted for convenient and 
reliable fits of weak bursts. 

The uncertainties in the interpolation of the background is an 
additional source of systematic errors which can be substantial 
for weak GRBs. This bias is difficult to estimate. It can be re- 
sponsible for large (more than \a) deviations for weak events. 
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Fig. 14. — Fitting count rate time profiles for the eleven weakest GRBs 
detected in our scan. The burst 10330a can hardly be recognized on this scale 
due to its long duration (~ 150 s) and smoothness. Nevertheless, this is a con- 
fident event clearly visible on a longer time scale and at a wider time binning. 
The time profiles have been shifted vertically for clarity. 

TABLE 1 
Parameters for the 1 1 weakest GRBs in the sample 



ID 


c 


/ 


Sign. 


Si 
(deg) 


^4 
(deg) 


08730a 


0.049 


0.072 


17 


18 


37 


08913c 


0.049 


0.088 


11 


23 


43 


09201c 


0.045 


0.076 


15 


33 


70 


09486b 


0.045 


0.076 


14 


24 


47 


09882b 


0.047 


0.077 


9 


33 


61 


10237b 


0.047 


0.081 


10 


39 


84 


10330a 


0.045 


0.072 


8.5 


37 


67 


10353d 


0.038 


0.066 


16 


20 


41 


10757a 


0.047 


0.081 


7.5 


24 


75 


10798a 


0.044 


0.075 


9.5 


39 


84 


10856d 


0.038 


0.088 


14 


44 


60 



ID is the event identifier consisting of TJD plus an identifying letter; c is 
peak count rate in units of counts s"' cm"^; / is peak photon flux in units of 
ph s"' cm"^ ; significance is estimated using the residual x^ relative to the lin- 
ear background in units of cr = \/2N, where N is the number of the degrees of 
freedom; S[ and S4 are the sizes of the Icr and 4a confidence areas (see §4). 



6.2. Bursts near the Detection Threshold 

The weakest events of the sample have peak fluxes below 
0.1 phs"' cm"^. They have been detected due to their long du- 
ration using the 4 or 8 second trigger time scales. They are 
not necessarily the least significant events. One can get a feel- 
ing of the weakest part of our sample from Figure 14 and Ta- 
ble 1, where we present information on the 11 weakest events 
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detected in the scan and classified as GRBs. The events in Fig- 
ure 14 and Table 1 are marked by their names consisting of 
the TJD when the event was detected and an identifying letter. 
All are more than 20 s long (shorter events are not significant 
at these brightnesses), are significant, and have good x^ maps. 
The detection efficiency in this brightness range is ^ 0.05. 

6.3. The Data Archive 

The electronic archive** (Stern & Tikhomirova 1999) contains 
the following data: 

1 . Table with event identifiers, arrival times, coordinates, 
peak fluxes, and duration estimates (Tgo and A^5o, i.e., 
the number of 1 .024 s bins in original DISCLA data, 
where the signal exceeds 50% of the peak value). 
BATSE trigger numbers are given for triggered events. 
Events identified with the GRBs of the K98 catalog are 
marked. Events in proximity of the data gaps are also 
marked. See Table 2. 

2. Fragments of raw DISCLA data in FITS format 
covering the events. The typical time interval is 300 
seconds for shorter (< 100 s) events. For longer events 
the time interval is extended. 

3. Time profiles as the sum of the count rates in the two 
brightest detectors and in energy channels #2 and #3 
(the 50-300 keV band) in ASCII format. This is the 
simplest and the most compact representation for the 
time profiles. 

4. Time profiles of GRBs as the fitting count rate, Ca, 
with subtracted background, see §3 and §5.1, for the 
four LAD energy channels in ASCII format. This 
representation has a better signal/noise ratio than the 
sum of the count rates in any combination of detectors. 

5. All distributions presented in this paper in both 
graphical and numerical form. 

7. TESTS FOR NON-GRB CONTAMINATIONS 

All our tests for non-GRB contaminations were performed 
on the sample of 3678 events detected up to TJD 11499 (8.6 
years of observations). 

7.1. Possible Types of the Contamination 

We can divide the possible non-GRB contaminations of the 
sample into the following categories: 

1. Poisson or non-Poisson background fluctuations. The 
least significant GRBs of our sample are too significant to 
be reproduced by pure Poisson fluctuations. There are how- 
ever significant variations of the background due to ionospheric 
phenomena, occultations of astrophysical sources, and noise- 
generating sources. Such events do not pass either the Sa,- 
criterion or the "isolation" criterion (see §4). Nevertheless, we 
should admit the existence of some cases where a fluctuation 
may pass all tests. This kind of fluctuations is almost time- 
symmetric. Therefore we can test for such a kind of contami- 
nation by searching for "negative bursts" in sign-inverted count 
rate data. We performed such tests for 200 days of data records 
taken in different intervals over a few years representing vari- 
ous background conditions. The rate of triggers was practically 

^The archive is available at: littp://www.astro. su.se/groups/liead/grb_arcliive.litml 



the same as in the normal scan. However, only 2 events passed 
the objective criteria described in §4 (mainly they were rejected 
by the Sa criterion). Both events are within the "2cr area" of 
Cyg X-1 and both are not enough isolated, so they should be 
rejected subjectively. Thus we constrain the contamination of 
sign-symmetric background fluctuations to be less than ^ 0.5% 
of the events. 

2. Ionospheric events such as particle precipitation. Particle 
precipitation is a frequent phenomenon taking place at high lat- 
itudes: the satellite flies through a cloud of energetic particles 
emitting bremsstrahlung. Then a more or less time-symmetric 
hump of counts is recorded in all or most of the detectors. Such 
events are easily recognizable by comparable signals in detec- 
tors facing opposite directions and, correspondingly, by a bad 
X^ map (many local minima of x^ over a wide area of the sky) 
and a large residual x^ (above 0.1 of its initial value). If the 
precipitation occurs at some distance from the satellite path, 
the time profiles in different detectors do not match and again 
we obtain a bad location fit. The most dangerous events are 
distant ionospheric flares with fast (shorter than hundred sec- 
onds) variability. We observed such events recognizing them 
by "context", i.e., they occur at a high general ionospheric ac- 
tivity at high latitudes and appear in long series. 

3. Solar flares and known variable X-ray sources (Cyg X-1, 
bright X-ray transients, and known X-ray pulsars). Solar flares 
are not a considerable source of the contamination as they are 
much softer than GRBs. There is some small overlap in hard- 
ness. For this reason, we reject some events which could be soft 
GRBs if they are consistent in location with the Sun. 

4. Unknown astrophysical sources emitting bursts similar to 
GRBs. These could be Galactic accreting black holes, hard 
bursts of other X-ray pulsars, or unknown phenomena (i.e., 
huge stellar flares). We cannot discriminate such events indi- 
vidually if they overlap with GRBs in their characteristics. We 
can only sense such contaminations statistically. Such a statis- 
tical analysis is presented below. 



7.2. Hardness Ratios 

GRBs are on average harder than any other kind of events 
that can mimic a GRB. We plot several examples of non-GRB 
events of known origin together with GRBs in a hardness - 
brightness diagram. Figure 15, where the hardness ratio is that 
between the counts in the 50-300 keV and the 20-50 keV en- 
ergy bands. Outbursts of Cyg X-1 have an almost constant 
hardness ratio. It is slightly smaller than the average ratio for 
GRBs. Ionospheric flares have a wider dispersion and partly 
overlap with GRBs in their hardness ratios, while being softer 
on average. Bursts of X-ray pulsars are on average consider- 
ably softer than GRBs. However, their hardness distribution 
overlaps with that of GRBs and some hard well-isolated pulses 
could be included into the sample as GRBs. Large (tens of 
events) contamination of this kind would be revealed in the an- 
gular distributions of GRBs. The softest events we deal with 
are solar flares. Their overlap with GRBs in hardness is very 
small as seen in Figure 15. 

If the weakest part of our sample were considerably contam- 
inated with some non-GRB events, we could observe this in 
the hardness distribution of weak events (this is hardly so for 
Cyg X-1 which has almost the same hardness ratio as weak 
GRBs). However, the hardness ratios for the weakest events 
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smoothly extend the hardness versus brightness correlation of 
brighter GRBs as seen in Figures 15 and 16. Thus the statistics 
of the hardness ratios for the weakest bursts does not demon- 
strate any indications of a considerable contamination. 
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FlG. 15. — Hardness ratio (the integral of the counts in channels #2 and 
#3, 50-300 keV, to that in channel 1, 20-50 keV) versus the 50-300 keV peak 
flux for GRBs (dots) and for different kinds of non-GRB events found in the 
BATSE records. Crosses: Cyg X-1; squares: ionospheric phenomena; stars: 
X-ray pulsars (two objects), circles: solar flares. 
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Fig. 16. — Hardness versus the 50-300 keV peak flux for the GRBs de- 
tected in the scan. The hardness ratio is here the ratio of the time-integrated 
peak count rate in the BATSE LAD channel #3 (100-300 keV) to that in chan- 
nel #2 (50-100 keV). Solid curve shows the median hardness ratio for trig- 
gered GRBs, while the dashed curve shows the same for non-triggered GRBs. 
Circles represent the 1838 non-triggered GRBs, while dots represent the 2068 
BATSE-triggered GRBs detected in the scan. Weak triggered GRBs are harder 
than non-triggered ones in the same peak flux range because of short hard 
events. They are weak at 1024 ms resolution but are easily triggered at 64 
ms resolution. The general trends approximately agree with the results of Ne- 
miroffetal. (1994). 



7.3. Angular Distributions 

Each kind of non-GRB events has its specific non-isotropic 
angular distribution which can affect the angular distribution 
of events of the sample if the latter is contaminated with non- 
GRBs. Here we analyze several angular distributions of events 
in the sample. Each of them is sensitive to some kind of non- 
GRB contamination. Since the weakest part of the sample is 
more likely to be contaminated, we analyze some angular dis- 
tributions using the subsample of the weakest GRBs. The polar 
angular distributions for both real and test bursts are consistent 
with the BATSE sky exposure function (Paciesas et al. 1999) as 
seen in Figure 17. 
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Fig. 17. — Polar angular distributions of GRBs and test bursts. Smooth 
curves show the BATSE sky exposure function (Paciesas et al. 1999). The 
GRB declination is denoted by S. 
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Fig. 18. — Distribution of locations of GRBs relative to zenith. The arrow 
shows the mean position of the Earth's horizon. (The cosine of the horizon 
angle varied between about —0.31 and -0.38 due to the changing altitude of 
CGRO.) 
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The distribution of GRBs relative to the Earth's horizon is 
shown in Figure 18. The shape of the distribution is reason- 
able: an isotropic distribution above the horizon and a clear step 
at the horizon. There is no excess of events towards the horizon 
which could appear if the sample contains many misidentified 
events of terrestrial (ionospheric) origin. 

No excess of events towards the Sun has been found. Fig- 
ure 19 demonstrates a very wide separation between solar flares 
and GRBs in a plot of hardness ratio versus angular distance 
from Sun. One can conclude from Figure 19 that the contami- 
nation of the sample by solar flares is negligible. 
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Fig. 19. — Hardness ratio vs. angular separation from the Sun for GRBs 
(dots) and several solar flares (circles). 
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Fig. 20. — Angular distributions of detected weak (P < 0.2phs"' cm"^) 
bursts relative to Cyg X-1, the Galactic Center, and the Galactic pole. Thin 
line histograms are for weak test bursts, and thick line histograms - for weak 
real bursts. 



The distribution of weak GRBs and weak test bursts (P < 
0.2phs"' cm"^) relative to the direction to Cyg X-1 is shown in 
Figure 20. There is an evident deficit of both test and real bursts 
in the direction to Cyg X-1. This is natural as we missed many 
weak bursts in that direction when subtracting the signal from 
Cyg X-1. The depression is slightly deeper in the case of test 
bursts. The corresponding expected number of real bursts in 
the two last bins (a cone of half opening angle 26° around Cyg 
X-1) is 11, while the real number is 15. The difference is not 
significant, of course. The formal Icr confidence interval for 
contamination by Cyg X-1 is 4±5 events. Brighter outbursts 
of Cyg X-1 are very rare. 

Contaminating X-ray pulsars and variable X-ray sources 
would give an excess of weak bursts in the direction of the 
Galactic Center. Actually there is, on the contrary, a deficit 
of both real and test weak bursts in the direction of the Galactic 
Center A small fraction of this deficit is explained by the cov- 
erage function. A stronger effect appears as a result of worse 
background conditions in that direction, mainly due to Cyg X-1 
(^ 60° from the Galactic Center direction). 

The formal la confidence interval for a Galactic Center sub- 
population in our sample for a cone of 37° half opening angle 
in the direction of the Galactic Center (where both of the X- 
ray pulsars shown in Figure 15 are located) is -10 ± 10 events. 
Similarly we do not see any excess of weak events towards the 
galactic plane. 
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Fig. 21. — Normalized distributions of detected bursts and particle precipi- 
tation events over the latitude of CGRO at the time of detection. Upper panel: 
weak events (P < 0.2phs" cm ); lower panel: all events. Thick line his- 
tograms: real bursts; thin line histograms: test bursts; dotted line histogram: a 
sample of 100 particle precipitation events detected during different time inter- 
vals of observations. The asymmetric shape of the burst histograms is mainly 
due to the South Atlantic magnetic anomaly where any detection is impossible 
because of a very large ionospheric background. 

7.4. Dependence on the Satellite Latitude 

The intensity of ionospheric activity strongly depends on the 
satellite latitude. Particularly, the frequency of particle precip- 
itations sharply increases at high latitudes. Therefore the dis- 
tribution of the satellite latitude, when an event of our sample 
was recorded, is sensitive to the contamination of the sample 
by the ionospheric events. Then we would have an excess of 
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events detected at large Northen and Southern latitudes. The 
distribution of events along the satellite path should not be uni- 
form as the background conditions are substantially worse at 
a high latitude. As far as we have common conditions for the 
detection of real and test bursts, an excess of real events rel- 
ative to the test bursts can be estimated. Figure 21 shows the 
results of such a test for all and for the weakest events. There 
is no excess of real events relative to the test bursts at high lat- 
itudes. An approximate la upper limit for the contamination 
of the sample by events of ionospheric origin at latitudes above 
±20° (where more than 90% of the particle precipitation events 
occur, see Fig. 21) is 42 events (i.e., 1.1%) for the full sam- 
ple and 17 events taking into account only the weakest part 
of the sample where a misclassification is more probable (at 
P < 0.2phs-' cm-2, the 680 weakest GRBs). 

7.5. Summary of Contamination Tests 
Summarizing the results of the tests: 

1 . The contamination of the sample by solar flares is 
negligible. 

2. The Cyg X-1 contribution is estimated as 4 ± 5 events. 

3. The contribution of all sources concentrated within 37° 
around the Galactic Center direction is constrained to 
be -10 ±10 events. 

4. The least constrained source of contamination is 
produced by ionospheric phenomena because the 
corresponding angular distributions are the least deter- 
mined. However, even in this case the contamination is 
constrained to be at the level of ^ 1%. 

We should admit the existence of some non-GRB isotro- 
pic background of events of a different nature but looking like 
GRBs. Then this is rather an issue of the classification of GRBs, 
which is beyond the scope of this paper. 



8. THE logW-log/' DISTRIBUTION 

. 1 . The Efficiency Function and the Corrected logN- 
Distribution 



logP 



The efficiency function of the search is defined as the ratio of 
the number of detected test bursts to the number of test bursts 
applied to the data versus the expected peak count rate, c^. It 
is shown by the histogram in Figure 22. For many applica- 
tions, it is useful to have a simple analytical representation of 
the efficiency function. We find that this function can be well 
approximated by 



£(q) = 0.70<^ 1-exp 



(5) 



where two parameters Ce.o = 0.097 countss"' cm"^ and v = 2.34 
(see the smooth curve in Fig. 22) 

The result of a direct correction of the measured peak count 
rate distribution (Fig. 8 and the dotted histogram in Fig. 23) us- 
ing the efficiency function is shown in Figure 23 by the solid 
histogram. Such a correction is valid for reconstructing the 
logA^-logP distribution only at peak brightnesses that suffi- 
ciently exceed the threshold such that the measured and the ex- 
pected count rates are almost identical. There, the errors in the 
peak count rates are moderate and symmetric (see Fig. 6) and 



the efficiency changes smoothly. One sees that the corrected 
logA^-logP distribution extends down to lower brightnesses al- 
most straight without any indication of a turn-over. This result 
differs from previous interpretations of the data, which implied 
that the logA^-logP distribution smoothly bends down at low 
brightnesses (see the data points of KOO in Fig. 23). 
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Fig. 22. — The efficiency of the off-line burst detection defined as the 
fraction of the test bursts detected in our scan versus the peak count rate, P. 
The efficiency is normalized to the total number of bursts occurring above the 
Earth's horizon (all test bursts were sampled above the horizon). The smooth 
curve shows the fit given by equation (5). 
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Fig. 23. — The differential log W — log P distribution corrected for the ef- 
ficiency function in absolute units for all 3906 GRBs detected in our scan 
(solid histogram). The original uncorrected distribution is shown by dotted 
histogram. The corrected distribution for the 3300 events in our sample with 
a duration longer than 1 s is shown by crosses. The corresponding distribu- 
tion for the 2265 GRBs (of any duration) found by KOO is shown by circles 
(the KOO data were transformed to count rates by applying the factor 0.75, see 
§5.2 and Figure 7). The leftmost dashed data point is probably affected by the 
threshold bias arising from errors in the peak count rates. 
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Our whole sample is not homogeneous because it contains 
short events where the peak flux estimate in 1 s time resolu- 
tion is wrong. To eliminate the corresponding bias we excluded 
from the distribution all events consisting of one bin (i.e., where 
only one 1.024 s bin is above 0.5 of the peak flux value). The 
corresponding efficiency function was calculated using only 
those test bursts that produce events satisfying the above cri- 
terion. Then we have a more or less homogeneous estimate of 
the peak count rates. The logA^-logP distribution of only long 
events is shown by crosses in Figure 23. Such a distribution is 
more relevant for cosmological fits than the total distribution. 

8.2. The Efficiency Matrix and \ogN —\ogP Close to the 
Detection Threshold 

When the burst peak flux is close to the detection thresh- 
old, the difference between the expected and the measured peak 
count rates can be significant. Therefore, in order to reconstruct 
the log A^- log P distribution at low fluxes, we have to take into 
account the matrix of the probabilities of detecting a burst of a 
given expected count rate with a given measured count rate. Us- 
ing test bursts from our sample, we constructed such a matrix. 
The "low brightness" end of this matrix is presented in Table 3 
(see also Fig. 6). We can see the asymmetry of the matrix near 
the threshold (i.e., the off-diagonal elements are larger than the 
diagonal ones) which results from the selection bias discussed 
in §5.2. 

In spite of the fact that we had 11 100 test bursts, the statistics 
at the low brightness end of the efficiency matrix is relatively 
poor and the error there is large. We can improve the situa- 
tion by approximating the matrix by a smooth function of the 
expected and the measured count rates, Cg and Cm, and by ex- 
trapolating this function to smaller count rates. We approximate 
the efficiency matrix with a factorized expression where the ef- 
ficiency function, the peak brightness error, and the selection 
bias were fitted independently: 



F{Ce,Cm) = E{Ce)- 



1 



:exp 



log \Cffi/ CifiQj 
2^2 



(6) 



where E(ce) is given by equation (5), the log-normal factor 
describes the relative error of the measured count rate, a = 
0.09(0. 08 /ce)'''^, and the selection bias is crudely expressed as 
Cm,o = Ce + 0.05 exp(-Q/0.05). 

As was mentioned above, the correction using the efficiency 
function is not exact especially near the threshold. One should 
deconvolve the observed peak brightness distribution (Fig. 8) 
using the efficiency matrix. This is easier to do with a forward 
folding method, i.e., fitting the observed distribution using a 
convolution of a hypothetical log A^- log P distribution with the 
efficiency matrix. We cannot use "data points" like those in 
Figure 23 for the hypothetical logA^-logP distribution as the 



forward folding in this case allows large fluctuations between 
neighboring data points and the result is unstable. One should 
use some smooth function as a hypothesis for the logA^-logP 
distribution, convolve it with the efficiency matrix (6) and com- 
pare with the data. Such a fit is beyond the scope of the present 
paper. 

9. CONCLUSIONS 

With a careful off-line scan of the daily archival data, BATSE 
becomes a more sensitive instrument with better known char- 
acteristics. The detection threshold changes by a factor 2, from 
^ 0.2phs"' cm"2 to ^ 0.1 phs"' cm"^. The total number of de- 
tected GRBs increases by 70 per cent. If we consider only long 
events (short GRBs could be a different phenomenon), the gain 
is a factor of ^ 2. The recognition of the weakest GRBs is still 
confident and the contamination with events originating from 
trivial kinds of non-GRBs is small. 

As a first result, using the new sample of GRBs with known 
detection efficiency, we pointed out that the log A^- log P dis- 
tribution does not show a turn-over at lower brightnesses im- 
plied by previous studies. One simple consequence is that 
the estimate of the number of GRBs in the visible Universe 
should be increased. Just the "visible" part of the logA^-logP 
distribution implies 1200-1300 GRBs occurring per year at 
peak fluxes exceeding 0.1 phs"' cm"^ (previous versions of the 
logA^-logP distribution implied ^ 600 GRBs per year above 
0.1 8 phs"' cm"2). A possible extrapolation of the new logA^- 
logP distribution to lower brightnesses would probably imply a 
much larger rate of up to several thousands of GRBs per year. 

The best possible efficiency and scan quality have still not yet 
been achieved. When we repeated the scan for a fraction of data 
scanned early on, we found 13 new GRBs per 60 days. This 
means that an additional considerable increase of the statistics 
of the weakest GRBs is possible (as is an extension of the statis- 
tics of the weakest test bursts). There are obvious possibilities 
for improving the test burst method as was discussed in §2 and 
in Stern et al. (2000b), for improving the approximation of the 
background, etc. This means that a new scan is desirable. 
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TABLE 2 
GRB CHARACTERISTICS 



ID" 


Time" 


BATSE Trigger'^ 


Peak Flux" 


R.A.'' 


Decl.= 


<5i = 


Ton' 


^50*^ Kommers 


Gap 




(s) 


or 1 if not 


(phs"'cm"2) 


(deg) 


(deg) 


(deg) 


(s) 


(bin) Catalog^ 


Flag' 


08369d 


71474 


107 


0.373 


220.0 


10.5 


6.1 


207 


37 





08370b 


6141 


1 


0.248 


68.4 


-6.4 


8.6 


83 


21 





08370f 


71005 


108 


0.178 


209.8 


-45.6 


34.1 


4 


2 





08371a 


2208 


109 


3.725 


89.0 


-16.0 


0.2 


95 


27 


1 


08371c 


20208 


110 


0.454 


341.8 


25.0 


7.9 


549 


6 





08372a 


5386 


1 


0.336 


246.2 


65.4 


7.6 


8 


5 





08372d 


80024 


111 


0.502 


78.1 


-22.5 


3.8 


84 


25 





08373a 


14897 


1 


0.257 


32.4 


-46.1 


21.4 


85 


5 





08373c 


32689 


114 


0.540 


76.5 


-19.2 


11.1 


370 


5 





08375b 


11481 


121 


1.324 


173.2 


9.1 


1.9 


87 


13 





08375c 


21537 


1 


0.259 


315.3 


-48.9 


17.2 


50 


13 





08376d 


61501 


130 


3.561 


133.6 


1.9 


1.4 


203 


15 





08377a 


30013 


133 


0.621 


124.3 


-1.9 


6.1 


174 


22 





08378b 


40667 


138 


0.316 


126.0 


-21.4 


17.1 


3 


1 





08379b 


25453 


143 


37.809 


88.9 


35.8 


0.0 


56 


3 





08379c. . . 


26090 


1 


0.291 


214.7 


-19.6 


8.2 


101 


13 






Table 2 is published in its entirety in the electronic edition of Tlte Astrophysical Journal. See also ftp://ftp.astro.su.se/pub/head/grb/catalogs/etable2.txt. 

A portion is shown here for guidance regarding its form and content. 

"Event identifier consisting of TJD plus an identifying letter. 

''Start of the event within TJD. 

"^BATSE trigger number or 1 for non-triggered events. 

''Peak flux in the 50 - 300 keV range. 

■^The best fit coordinates (J2000) and the size of the la confidence area. 

'Duration. 

^Number of 1.024 s bins in the original DISCLA data where the signal exceeds 50% of the peak value. 

''Flag identifying overlap with the K98 catalog. A 1 indicates existence in K98. 

'Flag characterizing the data quality with non-zero values indicating proximity to the data gap. 



TABLE 3 
The "low brightness" part of the efficiency matrix 



logc, 








logCc 










-1.3 


-1.2 


-1.1 


-1.0 


-0.9 


-0.8 


-0.7 


-0.7 


0.0 


0.0 


0.0 


0.0 


0.7 


11.7 


43.1 


-0.8 


0.0 


0.0 


0.3 


1.0 


12.1 


36.5 


11.0 


-0.9 


0.0 


0.0 


1.8 


10.4 


28.6 


11.2 


0.7 


-1.0 


0.0 


0.8 


6.8 


15.9 


9.5 


1.3 


0.4 


-1.1 


0.0 


4.3 


6.9 


6.5 


1.0 


0.1 


0.0 


-1.2 


0.0 


0.4 


2.9 


0.8 


0.2 


0.1 


0.1 


-1.3 


0.0 


1.2 


0.0 


0.0 


0.0 


0.0 


0.0 



Probability F(Ce,Cm) (in %) for the detection of a GRB with an expected peak count, c^, and a measured peak count rate, 
horizontal bins are uniform in the logarithm of c. The lower edges of the bins are given in terms of log(Cf ,„,). 



( counts s cm ). The vertical and 



